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We study the competing phases and the phase transition phenomena in an effective spin- 
orbital coupled model derived for pyroxene titanium oxides ATiSi2 0e (^4=Na, Li). Using the 
mean-field-type analysis and the numerical quantum transfer matrix method, we show that the 
model exhibits two different ordered states, the spin-dimer and orbital-ferro state and the spin- 
ferro and orbital-antiferro state. The transition between two phases is driven by the relative 
strength of the Hund's-rule coupling to the onsite Coulomb repulsion and/or by the external 
magnetic field. The ground-state phase diagram is determined. There is a keen competition 
between orbital and spin degrees of freedom in the multicritical regime, which causes large 
fluctuations and significantly affects finite-temperature properties in the paramagnetic phase. 
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1. Introduction 

It has been recognized that the orbital degree of free- 
dom plays an important role in transition metal com- 
pounds. 1-3 The interplay among the orbital, spin, and 
charge degrees of freedom produces a rich variety of phe- 
nomena. In particular, quantum and thermal fluctuations 
in the competition between the spin and orbital exchange 
interactions have attracted much interest. 

Pyroxene titanium oxides ATiSi2C>6 (j4=Na, Li) are 
typical examples of the systems where the interplay be- 
tween orbital and spin is expected. 4 In these insulating 
materials, each magnetic Ti 3+ cation has one d electron 
in t2 g levels. The lattice structure consists of character- 
istic one-dimensional (ID) chains of skew edge-sharing 
TiC>6 octahedra, which are bridged and well separated 
by SiC>4 tetrahedra. Hence, we can regard the system as 
a quasi ID spin- 1/2 system. 

These pyroxene titanium oxides exhibit a phase transi- 
tion to the low-temperature phase with a spin gap and a 
lattice dimerization at T ~ 200K, 4 ' 5 which reminds us of 
a spin-Peierls transition. However, the magnetic suscep- 
tibility shows a peculiar temperature dependence clearly 
distinct from that of usual spin-Peierls compounds. The 
magnetoelastic effect is not sufficient to explain the data. 
Therefore, the other scenario is needed to capture the 
physics of the transition. 

In the previous study, we derived an effective spin- 
orbital coupled model for the pyroxene compounds in 
the strong correlation limit of a multiorbital Hubbard 
model, and investigated it using the numerical method as 
well as mean-field (MF) analyses. 6 We have shown that 
the interplay between orbital and spin degrees of free- 
dom plays a decisive role in the spin-singlet formation. 
In particular, we have found that the system exhibits a 
nontrivial feedback effect between the orbital and spin 
correlations: At high temperatures, both correlations are 
antifcrro (AF) type and compete with each other. When 
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the Hund's-rule coupling is small enough compared to the 
onsite Coulomb repulsion, the AF spin correlation wins 
the competition and grows rapidly as temperature de- 
creases. The enhancement of the spin correlation changes 
the orbital correlation from AF to ferro (F) type, and 
eventually the F-type orbital correlation causes the spin- 
singlet formation. The experimental data are explained 
by the effective model semiquantitatively. 

Besides these results, we have also found that there 
is a chance to realize another type of ordered state in 
the effective model. Namely, if the Hund's-rule coupling 
is relatively large, the AF orbital correlation may over- 
come the AF spin correlation, yielding the spin-F ordered 
state with AF-type orbital ordering at low temperatures. 
This result indicates that the model may exhibit a phase 
transition between the spin-singlet and orbital-F phase 
and the spin-F and orbital-AF phase by changing the 
strength of the Hund's-rule coupling. 

The aim of this paper is to study the phase transition 
in the effective spin-orbital model in detail. Using the 
MF analysis and the numerical quantum transfer matrix 
(QTM) method, we clarify the nature of the ground-state 
phase transition and the spin-orbital interplay at finite 
temperatures in the critical regime. Another purpose of 
this study is to explore the effect of an external mag- 
netic field on the spin-orbital interplay. We find that the 
magnetic field also induces the same-type transition. 

The paper is organized as follows. In Sec. 2, we in- 
troduce the model Hamiltonian. The MF results for the 
ground state are shown in Sec. 3. There, we also discuss 
a basic mechanism of the phase transitions. In Sec. 4, 
we present the numerical results at finite temperatures 
obtained by the QTM method. Section 5 is devoted to 
summary and concluding remarks. 

2. Model Hamiltonian 

The microscopic model for the pyroxene compounds 
has been proposed in Ref. 6. Here, we briefly overview 
the derivation. 
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Fig. 1. Schematic picture of the ID multiorbital Hubbard model 
for the pyroxene titanium oxides. The ellipses represent Ti sites 
with doubly-degenerate orbitals, and the gray curves show the a 
bonds. 
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We start from a ti g multiorbital Hubbard model. In 
the pyroxene compounds, due to the skew chain struc- 
ture, the TiOe octahedra are distorted and the threefold 
tig levels are split into a low-lying doublet and a sin- 
gle higher level. 5 It is therefore reasonable to consider 
the ID Hubbard model with doubly-degenerate orbitals. 
For the hopping term, considering the edge-sharing net- 
work of TiOe octahedra, we take account of only the 
cr-bond transfer integrals, i.e., the transfers between the 
nearest-neighbor (NN) pairs with the same orbitals lying 
in the same plane, and neglect other small matrix ele- 
ments. Due to the peculiar zigzag structure of the TiOe 
chain, 4 ' 6 the NN transfer integrals i"f +1 , where i is site 
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where n is an integer. The schematic picture of the a 
bonds is shown in Fig. 1. As for the onsite interactions, 
we consider the intra- and inter-orbital Coulomb repul- 
sions U and U' , and the Hund's-rule coupling Jh, which 
satisfy the relation U = U' + 2Jh- 

From the second order perturbation in the strong- 
coupling limit to- <C U, U' , Jh, the effective spin-orbital 
Hamiltonian (Kugel-Khomskii type 1 ) is obtained as 



h oAF 



a i,i+l 
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where Si is the S — 1/2 spin operator while Tj = ±1/2 
is the Ising isospin which describes two orbital states at 
site i. Note that the orbital isospin interactions are of 
Ising type because the transfer integrals (1) are orbital 
diagonal and only one of t}^ +1 or tf 2 i+1 is nonzero on 
each bond. In Eq. (2), we included the term for the Zee- 
man coupling with the external magnetic field h for later 
use. The coupling constants in Eqs. (2)-(4) are given by 
parameters in the original multiorbital Hubbard model 
as J = {t a ) 2 /U, A = 3/4(1 - 377) + 1/4(1 — 77), B = 
1/(1 -37?)- 1/(1 -v), and C = 4[l/(l + r ? ) + 2/(l-y 7 )]/3 
where rj = J-r/U . Hereafter, we set J = 1 as an energy 



Fig. 2. Schematic picture of (a) the sD-oF state and (b) the sF- 
oAF state. The rectangles in (a) represent the spin-singlet pairs. 
The ground-state energies per site are also shown. 



unit and use the convention of the Boltzmann constant 
fc B = 1. 

We note that the JT-type orbital-lattice couplings are 
also present in real compounds. For the pyroxene com- 
pound NaTiSi20g, the JT stabilization energy is esti- 
mated to be ~ 90K, which is considerably smaller than 
J ~ 200-300K. 6 In this paper, for simplicity, we consider 
the ID spin-orbital Hamiltonian Tt so only. Indeed, for the 
realistic value of the JT coupling, the low-temperature 
properties are dominated by the spin-orbital interplay 
intrinsic to Ti so , and the JT coupling just assists to es- 
tablish the long-range orderings at finite temperatures. 6 

3. Mean- field analyses of ground-state phase 
transitions 

3.1 At zero magnetic field 

In this section, we consider the ground-state proper- 
ties of the model (2) at zero magnetic field h — by 
using a MF-type argument. As discussed in the previ- 
ous study, 6 we may expect two different types of long- 
range ordered states in the realistic range of parameters. 
When the orbital ordering is F-type, e.g., (Ti) = 1/2 
for all i, Eqs. (3) and (4) show that the chain is dis- 
connected at every other bonds. The spin-exchange in- 
teractions for remaining isolated NN pairs are antifer- 
romagnetic and form the two-site spin-singlet states, re- 
sulting in the spin-dimer and orbital-F (sD-oF) ground 
state [see Fig. 2(a)]. On the other hand, if the orbital 
ordering is AF type, e.g., (Ti) = (— l) l /2, the system 
becomes a uniform ferromagnetic spin chain. Hence, the 
ground state is the spin-F and orbital- AF (sF-oAF) state 
[Fig. 2(b)]. Comparing the ground-state energies shown 
in Fig. 2 at h = 0, wc find that the ground state is sD-oF 
for 77 < while sF-oAF for i] > r/ c . The critical point 
is ?7 C = (v73 — 8)/3 ~ 0.181. This transition controlled 
by 77 is of first order and accompanied with finite jumps 
of the spin gap (from JC to 0) and the magnetization 
(from to the saturated value ±1/2). 

It should be noticed that the ground-state energies ob- 
tained above fully include quantum fluctuations since the 
orbital interactions are classical Ising type and since the 
spin ordering is fully polarized F type (no quantum fluc- 
tuation) or the two-site singlet state (whose energy is 
— JC per singlet pair exactly). In this sense, these re- 
sults are exact unless an ordering with a period longer 
than two-site spacing takes place. We will comment on 
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Fig. 3. Schematic picture of typical processes of the virtual hop- 
pings in (a) the sD-oF state and (b) the sF-oAF state. The energy 
gains for each process are also shown. 



the possibility of the ordering of longer period in Sec. 5. 

3.2 At finite magnetic fields 

The transition between the sD-oF and sF-oAF phases 
is also driven by the external magnetic field h. The mag- 
netic field suppresses AF spin correlations and favors 
spin-F state, and hence we may expect a field-induced 
transition from the sD-oF phase to the sF-oAF phase. 
Comparing the energies in Fig. 2 including the Zeeman 
contribution, the critical field h c is calculated as 

2(3 - 16t? - 3r? 2 ) 
c 3(1- 3^(1 + V)(l-V)' [> 
This is also the first-order transition; at the critical point, 
all the isolated singlet pairs collapse at once, leading to 
the fully polarized state. We note that the result (5) is 
also exact in the same sense as that in Sec. 3.1. The phase 
diagram is shown in Fig. 6 in comparison with numerical 
results. 

3.3 Competition between AF superexchange interaction 
and Hund's-rule coupling 

The phase transitions discussed in the preceding sec- 
tions can be viewed as a result of the competition 
between the spin-superexchange interaction and the 
Hund's-rule coupling. The following discussion is sketchy 
but useful to illustrate the competition. For simplicity, 
we consider the case without the magnetic field although 
the discussion is applicable to the case with a finite h. 

Figure 3 shows typical contributions in the perturba- 
tion processes to derive the model (2): (a) shows a typical 
process for the sD-oF configuration, which contributes 
to the AF spin-superexchange interaction. The energy of 
the intermediate state is U. On the other hand, (b) shows 
a typical process for the sF-oAF configuration, whose in- 
termediate energy is U' — Jh = U — 3Jh- Hence, the 
Hund's-rule coupling Jh > tends to enhance the energy 
gain of the latter process and stabilize the sF-oAF state. 
We note that the sF-oAF state is indeed favored for any 
positive Hund's-rule coupling (77 = Jh/U > 0) in the 
two-orbital model with uniform transfer integrals, i.e., 
r;ii +1 = /j? 2 i+1 = t a > for all NN bonds. 7 - 8 On the con- 
trary, in the present system, the energy gain by the su- 
perexchange process (a) is enhanced in the sD-oF state: 
Since the alternating transfer integrals (1) break up the 
system into independent two-site spin-singlet pairs, the 
ground-state energy for each connected bond is amplified 
up to — JC. Moreover, in the AF orbital configurations 
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Fig. 4. (a) Nearest-neighbor spin correlation Cs and (b) the 
nearest-neighbor orbital correlation Co for 0.160 < rj < 0.190. 



(b), there is no energy gain from the quantum singlet na- 
ture because the orbital isospins are of Ising type. These 
facts give rise to a finite regime of the sD-oF phase for 
the small values of r\. 

4. Numerical results at finite temperatures 

In this section, we present our numerical results for 
thermodynamic properties of the model (2) obtained by 
the QTM method. 9, 10 Details of the method is described 
in Ref. 6. In the following, we calculate the NN spin and 
orbital correlations, C s = (Si ■ S^+i) and C Q = (TiT i+ i), 
respectively, to clarify the spin-orbital interplay at fi- 
nite temperatures and the instabilities toward the ground 
states. The data are for the Trotter number M = 4 and 
we have checked the M-convergence of the data. 

4-1 At zero magnetic field 

First, we show the numerical results in the absence of 
the magnetic field. Figure 4 shows the data of C s and C 
for several values of 77. At high temperatures, both spin 
and orbital correlations are negative, i.e., AF type, and 
compete with each other. For 77 < 0.180, the AF spin cor- 
relation wins the competition and yields the sign change 
of the orbital correlation from AF to F type. After the 
sign change, both C s and C grow cooperatively as tem- 
perature decreases, and eventually appear to approach 
C s = —3/8 and C — 1/4, the values expected for the 
sD-oF ground state. On the other hand, for 77 > 0.182, 
the AF orbital correlation wins the competition and the 
spin correlation is changed from AF to F type. The cor- 
relations appear to converge to C s — 1 /4 and C Q = — 1/4 
as T — > 0, indicating that the system shows the sF-oAF 
ground state. These results suggest that the ground-state 
phase transition between the sD-oF and sF-oAF phases 
occurs at 77 = ?7 C ~ 0.181, being consistent with the result 
of the MF analysis. 

It is remarkable that for 77 close to rj c the competition 
between the orbital and spin correlations remains serious 
down to very low temperatures T ~ 0.1; both the correla- 
tions are gradually suppressed as temperature decreases, 
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Fig. 5. (a) Nearest-neighbor spin correlation C B and (b) the 
nearest-neighbor orbital correlation Co for r\ = 0.170 and sev- 
eral typical values of h. 
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Fig. 6. Ground-state phase diagram in the 77 versus h plane. The 
solid curve represents the mean-field phase boundary [Eq. (5)]. 
The circles correspond to the transition points estimated by the 
numerical calculations. 



and the convergence to the values expected for the or- 
dered states is not so clear. This peculiar behavior of the 
correlations indicates that the keen competition between 
orbital and spin severely affects the finite-temperature 
properties in the multicritical regime 77 ~ r) c . 

4-2 At finite magnetic fields: phase diagram 

In Fig. 5, we show the results of C s and C a for finite 
magnetic fields at r] = 0.17. Similarly to the results in 
Sec. 4.1, there occurs the transition from the sD-oF phase 
to the sF-oAF phase by increasing h, and near the critical 
field the severe spin-orbital competition survives down 
to very low temperatures (see the data for h = 0.3). We 
estimate the critical field as h c ~ 0.25 for r\ = 0.17, being 
consistent with the MF result, h c ~ 0.27. 

Performing the similar calculations for various values 
of 77, we obtain the ground-state phase diagram in the 
plane of 77 versus h in Fig. 6. The estimates of the crit- 
ical field by the QTM method agree well with the MF 
prediction. 



5. Summary and concluding remarks 

We have studied the effective spin-orbital model with 
alternating couplings [Eq. (2)], which is derived for the 
pyroxene titanium oxides ATiSi20g (A = Na, Li). Using 
the MF-type analysis and the numerical QTM method, 
we have shown that the model exhibits the transition 
from the spin-dimer and orbital-F phase to the spin-F 
and orbital-AF phase as the Hund's-rule coupling pa- 
rameter 77 = Jh_/U increases. The critical value of i] is 
estimated to be 77 c ~ 0.181. Moreover, we have also in- 
vestigated the effect of the magnetic field h on the phase 
transition and determined the phase diagram in the 77 
versus h plane. The results indicate that the field-induced 
phase transition occurs in the system with rj < n c . We 
have found a severe competition between orbital and spin 
degrees of freedom remaining down to very low temper- 
atures near the phase boundary. 

As mentioned in Sec. 3, the present MF analysis of the 
ground state is exact as far as we consider the orderings 
of two-site period. The QTM results, which can include 
longer-period orderings, appear to agree well with the 
MF estimates of the ground-state phase boundary. This 
agreement suggests that there is no intermediate phase. 
However, as we noted in Sec. 4, in the critical regime 
the QTM results are not so conclusive to determine the 
ground state because of the lack of very low-temperature 
data; there the QTM method suffers from systematic er- 
rors due to the finite Trotter number, and the accuracy is 
assured for T > 0.07 in the present study. Therefore, it is 
not ruled out to have novel phases in the narrow critical 
region. This interesting issue is left for future study. 

It is a challenging problem to observe the field-induced 
transition experimentally The critical field h c becomes 
large as 77 deviates from ?7 C ; from the estimate of J = 
200-300K, h c reaches 30-50T at 77 = 0.175. Therefore, to 
observe the transition in the experimentally accessible 
range of the magnetic field, it is needed to synthesize a 
compound with i] smaller than but close to rj c . Once an 
appropriate compound is prepared, one may observe an 
abrupt change of the magnetization at h c , from zero to 
the saturated value. 
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